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Recent lattice calculations have confirmed that QCD static potential for sources in dif- 
ferent representations of the gauge group is proportional to eigenvalue of the corresponding 
<—{ , quadratic Casimir operator with an accuracy of a few percents. We review the present the- 

oretical status of the "Casimir scaling" phenomenon and stress its importance for analysis 
£f~^ | of nonperturbative QCD vacuum models and other field theories. It is argued that Casimir 

scaling strongly advocates the property of Gaussian dominance and we propose different 
£NJ ' lattice tests to improve our understanding of these phenomena. 
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I. INTRODUCTION 



There is little doubt nowadays that QCD is the true theory of strong interactions and in perturbative 
domain agreement between theory and experiment is impressive. At the same time for nonperturbative (NP) 
phenomena an exact analytic formalism is still lacking and one is dealing with a set of models or approaches, 
which are not derived directly and rigorously from the QCD Lagrangian without any assumptions and most 
of these models are introduced ad hoc. 

Comparison of model predictions with experiment allows to select models describing one or another set of 
phenomena, e.g. constituent quark model works reasonably well for hadron spectrum, while parton model 
does the job for high-Q events and instanton model is suitable for describing chiral effects. However these 
models are not well connected to each other and may contradict other phenomena, e.g instanton gas model 
lacks confinement. In addition there is a class of specific models with a special purpose to explain vacuum 
structure of QCD with its property of confinement: magnetic monopoles with abclian dominance hypothesis, 
center vortices model etc. 

Somewhat apart from particular microscopic models a nonperturbative approach is being developed for 
the last ten years |Q (see also review Q and references therein), so called field correlator method. It is based 
on QCD and expresses all observables ultimately in terms of gauge-invariant field strength correlators. The 
crucial assumption of the method which makes it workable is that one can restrict oneself to a few lowest field 
correlators; taking only two-point Gaussian correlator one has Gaussian stochastic model of QCD vacuum. 
Extracting Gaussian correlator from lattice data ^ one is able to predict confinement property of QCD and 
calculate quantitatively an impressive amount of data on hadron spectra and other properties |^,^,^| • 

It is clear that QCD vacuum structures derived from different models should not necessarily be compatible 
with each other and also with GSM, while they can describe equally well some features of the vacuum 
dynamics, for example actual value of string tension. Taking gross features of the vacuum one should 
distinguish two different scenarios: one, based on an ensemble of classical solutions, i.e. of coherent lumps 
of fields may be called a coherent picture of the vacuum. Here belong instanton model, effective abelian 
models, central vortex models and some other models. On the opposite side there is Gaussian stochastic 
model where stochastic properties of the vacuum are taken to their extreme, so that all higher field correlators 
are suppressed and the picture of vacuum quantum correlations is close to that of the Gaussian white noise. 
It would be of great importance to understand first of all what the real QCD vacuum is, to which of these 
two different pictures it belongs. 
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The most important theoretical tool to treat NP QCD available up to now is numerical simulation of 
the theory on the lattice. Confinement and other NP phenomena have been extensively studied in this 
approach. Data obtained in this way can be compared with results of real experiments as well as with 
theoretical predictions, based on particular NP models. The latter comparison is of special interest. Since 
the only essential dynamical input of the lattice theory is QCD Lagrangian such comparison is important to 
distinguish between models, which are in agreement with QCD (at least in lattice formulation) and those, 
which are not. 

The property of confinement in QCD has two important facets. From observational point of view it exhibits 
itself as absence of particles, carrying nonzero color charge as asymptotic states in Nature (in particular, 
massless vectors and light spinors) . On the lattice there is another criteria of confinement - area law: there 
exists constant force (equal to » 15 tons in SI units) between static fundamental color charge and anticharge 
at large distances (if no dynamical charges can be created from the vacuum to screen it). Most of NP models 
which can be compared with lattice measurements are designed to describe confinement of colour charge and 
anticharge in the fundamental representation of the gauge group SU(3), i.e. the area law for the fundamental 
rectangular Wilson loop and hence linearly rising potential between static quark and antiquark. There is 
a lot of information about static fundamental potential obtained from lattice simulations both in quenched 
and unquenched approximation. Being of great value, these data cannot answer the questions posed above 
and one needs another type of lattice experiments which were done recently with great accuracy |6|-||] . These 
are the investigation of interaction between static charges in higher SU(3) representations. In this way one 
can extract information about field content of NP QCD vacuum, which is not available if only fundamental 
charges are considered. 

The studies in this direction have a long history (see |7j and references therein). Recently two sets 
of accurate lattice measurements of the static potentials for fundamental and higher representations of 
the gauge group SU(3) have been presented by G.Bali MM and independently by S.Deldar ||. The data 
strongly support the so called Casimir scaling hypothesis j^] which states the proportionality of the static NP 
potentials for different representations to eigenvalue of the quadratic Casimir operator for the corresponding 
representation. Actually, the level of violation of Casimir scaling (CS) for the string tension is found to 
be at the range 5 — 15% in the measurements H and at the level of a few percents in the results J(|, 
while extrapolation to the continuum limit performed in remarkably indicated the violation of CS law 
not exceeding statistical errors, of the order of 1%. All that give a hint that CS is indeed a fundamental 
property of QCD (at least on the lattice). Surprisingly, however, the analysis of Jl(|[ll| has shown that to 
incorporate such simple feature in some natural way is not an easy task for most popular QCD vacuum 
models. Let us make a comment about general reason for that. In most of coherent QCD vacuum models 
the corresponding field configurations are endowed with nontrivial colour structure, whose gauge-invariant 
meaning is often not transparent. The CS, on the other hand, means that at least for the particularly chosen 
NP quantity - static potential - only the simplest invariant of the given representation, eigenvalue of the 
quadratic Casimir operator, i.e. just square of the (nonabelian) charge characterizes the potential. Thus 
strongly interacting theory with complicated vacuum paradoxically demonstrates the same scaling law of 
the potential as it would be in perturbation theory at tree level. In terms of the confining field configuration 
ensemble it implies either delicate (and as we will try to advocate in what follows, rather unnatural) fine 
tuning of the contributions coming from more complicated color structures (higher powers of the quadratic 
Casimir and higher Casimirs) or actual smallness of such structures. The phenomenon of CS seen from 
different NP QCD models points of view will be discussed in Section III of the present paper. 

It should be said that many ideas connected with the meaning of CS and extensively discussed in the 
present paper were formulated already in the original papers || , |l^] , |L3| . Nevertheless we find it useful to 
collect them here all together with new developments since the results from [Q demonstrating much higher 
level of accuracy of the phenomenon are calling for reconsideration of the theoretical background behind it. 

The aim of the present paper is to draw possible consequences from the data on Casimir scaling constructing 
the most logical scheme of the QCD vacuum structure, respecting CS. The paper is organized in the following 
way: Section II is devoted to definitions and general discussion of the lattice data on CS, analysis of different 
field theoretical models from CS point of view is presented in Section III, Section IV contains our proposals 
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of additional observables aimed to clarify the CS behavior and structure of the QCD vacuum, and Section 
V presents conclusion and outlook. 



II. GENERAL FORMALISM 



The average of the Wilson loop W(C) for the rectangular contour C = R x T (which we choose in the 
'34" plane) for the dimension D representation of the gauge group SU(3) is given by 

(Wn(C)) = (ti d Pexp ( ig J Ap? a dz}\ ) (1) 

Here the normalized trace Tr d is defined as Tr £>1 = -^Tr 1 = 1, fundamental generators are normalized 
according to Tr T a T b = 5 ab /2. The SU(3) representations labelled by D = 3, 8, 6, 15a, 10, 27, ... are charac- 
terized by 3 2 — 1 — 8 hermitian generators T a which satisfy the commutation relations [T a ,T b ] = if abc T c . 
One of the main characteristics of the representation is an eigenvalue Cd of quadratic Casimir operator C D , 
which is defined according to 

C [ d ] = S ab T a T b = T a T a = C D -l (2) 

Our normalisation is such that Cd — N c for adjoint representation of SU(N C ), it differs by factor N c from 
normalisation of |lj] . Since any simple algebra of rank k has exactly k primitive Casimir-Racah operators 
p5| , p^| of order mi, ..,mfc, it is possible to express those of higher order in terms of the primitive ones. In 
the case of SU(3) the primitive Casimir operators are given by Cjp and — d a bcT a T b T c . The higher 
rank Casimir operators are defined as 

c« = rf|[;., r) r..T- (3) 

where the totally symmetric tensors df$ i ■, on SU (N > 3) are expressed in terms of dik and dijk (see, e.g. 
]14j and references therein). Following the notations from ^ we introduce Casimir ratio du — Cjj/Cp, where 
the fundamental Casimir Cp — (N^ ~ 1)/2N C equals to 4/3 for SU(3). Generally, for 5/7(3) representations 
labelled by Dynkin coordinates (fi, v) the eigenvalue of C D is to be calculated according to the following 
expression (see, e.g. @): 

C D (fJL,v) = i (/i 2 + + v 2 + 3/i + 3^) (4) 
o 

Notice the positivity of the quadratic Casimir charge (it is not true, for example, for the cubic one). 
The static potential between sources at the distance R in the representation D is defined as: 

V D (R) = - Hm hog (W(C)), (5) 

1 — >oo 1 

We are not considering the possibility of screening of the potential by dynamical degrees of freedom at the 
moment, see discussion below. The l.h.s. of ([|) can be decomposed into two pieces of different physical 
meaning: actual independent potential and self-energy part not depending on i?, which explicitly contains 
ultraviolet cutoff scale a. For the data extrapolated to the continuum limit (see |Q for details) one assumes 
the proper subtraction procedure to be applied and we will denote by Vjy(R) only the former part of the 
potential in this case. Notice, however that the numbers shown in Table 1 of the present paper are based on 
raw lattice data from Q where no any special subtraction have been performed (a comparison shows that 
with the discussed accuracy the self-energy part satisfies CS as well). 
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TABLE I. The Casimir-scaling and Casimir-violating string tensions and shifts ]Ll|. Based on the lattice data 
from H . All quantities with the hats are scaled according to u = u * 10 4 



D 




Aa< 4 > 








X 2 /dof 


8 


-3.5 


1.2 


-2.5 


2.8 


0.004 


19 / 43 


6 


-6.4 


1.2 


1.0 


2.6 


0.007 


26 / 42 


15a 


-5.2 


0.6 


-0.6 


1.1 


0.003 


39 / 42 


10 


-4.9 


0.5 


0.2 


1.0 


0.003 


22 / 41 



In the confinement phase lattice data are well described by the sum of perturbative Coulomb part, confining 
linear and constant terms: 

V D {R) = a D R + v D + ^§- (6) 

where all coefficients are D-dependent. The Coulomb term is now known up to two loops in the continuum 
limit |l^] and up to one loop on the lattice |ll| and is proportional to Cd in both cases. It is not trivial that 
perturbative interaction is proportional to Cd at two loops, and it is not clear what will happen with more 
loops taken into account. The CS hypothesis |9| states, that the NP confining potential is also proportional 
to the first power of the quadratic Casimir Cd, i.e. for the string tensions one should have (Jd/<Jf = do- 
As it was already discussed in the Introduction, this scaling law is in very good agreement with the results 
found in Earlier lattice calculations of static potential between sources in adjoint representation jl!| 

are in general agreement with [[|, however, deviations from scaling at the level of 10% are found in |2C|| , in 
particular, the value of is closer to 2 than to 9/4 in [^o) (see further references and discussion in H]). 

The lattice data can be understood in terms of bounds on higher Casimir terms in the potential (^) 10 
Namely, one can write quite general decomposition for the expression (g) of the following form: 

V D (R) = d D V^(R) + d 2 D V^(R) + ... (7) 

where the dots denote terms proportional to higher powers of do as well as to higher Casimir ratios. It 
is worth mentioning, that since there are two independent Casimirs in SU(3) case, the series in (^) is in 
fact Taylor expansion of a function depending on two variables (i.e. quadratic and qubic Casimirs) with R— 
dependent coefficients. Notice also, that for self-adjoint representations (those having (p,p) form in Dynkin 
coordinates) the cubic Casimir vanishes (see, e.g. |Q), so the potential depends only on do- 

The form of Taylor expansion (0) is strongly motivated by perturbation theory series as well as by cluster 
expansion (|2^). However it is not in one-to-one correspondence with any of them. For example, the lowest 
quadratic correlator {{FF)) contributes only to V^(R), while higher ones contribute to V^(R) as well as 
to higher terms. Neither the expression (R) has a meaning of perturbative series in do- 

The terms V (n) (i?) in the expansion (n) are representation independent. One can parametrize the first CS 
violating term in the same way as the total potential Vd(R) as V^(R) — + a^R, where we omit 
possible CS-violating contribution to the Coulomb potential since at small distances its magnitude should 
in principle be small due to asymptotic freedom and independent analysis confirms this. Here 
measure the ^^-contribution of the CS violating terms to the total potential (their artificial dependence on 
D could come from higher powers of dp, and should disappear if all terms in the expansion (Q) are taken 
into account). Notice, that we do not need to specify the coordinate dependence of V^ 2 \R). Some results 
based on the data from Q are presented in Table 1. The author of || used anisotropic lattice with the 
spatial unit a s — 0.082 Fm and anisotropy £ ~ 4. Standard x 2 fit was performed for the data in the whole 
range of all measured R, since no fingerprint of screening is seen up to the largest distances explored in |(| . 
Errors shown in Table 1 include statistical and systematical ones. 

Several comments are in order. First of all it is seen that the CS behaviour holds with a good accuracy, of 
the order of a percent with the reasonable x 2 /dof. The ratio /a^ is less than 1% for all cases presented 
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in Table 1. The value of the constant term t/ 4 ' is found to be compatible with zero within the error bars for 
all considered D, while it is not the case for <j( 4 K The value of for sextet, for example, is found to deviate 
from zero at the level of approximately five standard deviations. We have not found any strong systematic 
dependence of cr^ on D, which presumably confirms the validity of the expansion (Q) and shows, that the 
omitted higher terms do not change the picture in a crucial way. The negative sign of the correction 
is probably related in Euclidean metric with the fact, that the fourth order contribution is proportional to 
(ig) 4 > while the Gaussian term is multiplied by (ig) 2 < 0. 

Recently the continuum limit extrapolation of the lattice data on static potentials for different represen- 
tations was performed in 0. The actual ratio of the potentials Vd(R)/Vf(R) measured in demonstrated 
the deviation from CS (i.e. from a constant equal to do) at the level of statistical errors not exceeding a few 
percents. 

It is important to stress at the same time that CS property is not universal at all possible distances. The 
area law implies linear asymptotics of the potential at large R. There is however an important restriction: 
even in quenched approximation with no dynamical quarks included one should take into account the effects 
of string breaking. The conventional picture of this phenomenon suggests that in case of charges in higher 
representations the parts of the potential corresponding to the octet components of the sources are to be 
screened by dynamical gluons from the vacuum and so called gluelump states are formed eventually. In 
particular, zero triality (iV c -ality in general case) representations are screened completely. For example, in 
case of adjoint charges one has 8®8 = 27ffil0ffil0ffi8ffi8ffil with the singlet component. Screening for 
higher representations might require more dynamical gluons. In case of nonzero triality representations the 
noncompensated triplet is always present (for example, 6<g>8 = 3©6©15® 24) which makes asymptotic 
string tension to be equal to the fundamental one eventually. 

It is obvious that in the region where the string is broken the expressions of the form (^]), (^) have no 
sense. The modification of due to this effect was considered in |l3|,^l| in the strong coupling expansion, 
resulting in the following expression for adjoint Wilson loop in large N c approximation 

(W adj (C)) = exp[-a adj • Area(C)] + exp[-M ff ; • Perimeter(C)] (8) 

It follows from (||) that the potential (||) has two different regimes at R < R c and R> R c (R c is the critical 
distance, where the second term in (^) starts to dominate). Namely, one has in T — ► oo limit: 

V(R) = oR 6(R C -R) + oR c 6{R - R c ) 

The potential becomes flat after one reaches the critical distance. The estimates of the gluelump mass show 
p2| that the actual value of R c is rather large, of the order of 1.5Fm for the SU (3) gauge group in the adjoint 
case. In other words, the lightest gluelump is heavy in the units of a. At the distances actually explored in 
H H the second term in (||) contributes at the level of lattice statistical errors. This makes the "CS region" 
rather wide, at least for not very large D - from the smallest distances where tree level perturbation theory 
works well to the onset of screening. 

It is also worth mentioning that in the limit T — > oo the force is discontinuous according to (|^) at the 
point R c , switching from the constant value equal to a at R < R c to zero value at R > R c . This dynamical 
screening seems to be different from the center vortex kinematical screening scenario where the potential is 
becoming flat with the distance smothly (see below) even in the T — > oo limit. 

Let us show how (|§|) appears in dynamical background perturbation theory [f22f . We concentrate on the 
adjoint case for simplicity, generalization to other representations can be done in the same way. To this 
end we split gluon field as = + ga^ , where represent the confining background and a M - the 
valence gluon field. Referring the reader to [^|| for details, one gets the valence Green's function in the 
background Feynman gauge in the form G^(B) = (D 2 {B)8 lJiU + 2iF flv ) and the result of integration over 

valence gluons at the lowest order yields in the partition function the factor [DetG(B)] 2 . The averaging in 
(|), @ turns out to be 
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(W)i 



DetG(B)] * W(B) 



[DetG(B)]" 



(9) 



In a similar way quark loops are accounted for by the factor Det(l) + im) instead of [DetG[-B]]~ 5 in (|9|). 
The highest terms in ga^ expansion can be calculated systematically. 

The next step is the standard loop expansion of the determinant augmented by the world-line (Feynman- 
Schwinger) formalism : 



[DetG(B)]" 5 = exp 



OO 



= exp j {»'„,,, (£()} (10) 



where the curly brackets stand for the path integral over contours, forming the loop and the corresponding 
proper time integration. Expanding the exponent in ( |Io|) and keeping only the first two terms one has 



(W) = (W(B))j 



-((W(B){W adj (B)})) B 



(11) 



where double brackets ((..)) denote the connected correlator. The final step is the representation of the last 
term in (jll|) as a Green's function of two gluelumps, so one finally arrives to the asymptotic expression (|§|). 

Numerically, estimating M g i as 1.2 GeV for the adjoint source p2| , one gets for the critical distance 
R c 2M g i/a% as 1.2 Fm. Therefore as it was already stressed the preasymptotic region (where the confining 
potential is linearly rising for all D) is rather wide which allows to establish CS up to the distances where 
signal disappears into noise on the lattice. Additional measurements at larger distances could hopefully shed 
some light on the string breaking and physics of gluelumps and establish representation dependent bounds 
of the CS regions. 

Let us conclude this section with the remark about large N c limit. First of all, the "CS region" expands 
up to infinity in this case, as it is clearly seen from (||). It was also advocated in |2l| that the ratio d a dj is 
equal to 2 if the limit N c — > oo is taken (which is in line with both flux counting rule and CS). However we do 
not see clear arguments in continuum theory why CS should become exact if N c — ► oo in general case for all 
higher representations. It would be very interesting to establish in future lattice measurements the relative 
deviation from CS as a function of N c , at least for N c — 2. .5 (see also p3|). Notice, that the standard large 
N c factorisation arguments are not applicable to the cumulants entering (|22|) since they are color-irreducible. 
In particular, the proper understanding of Gaussian dominance property in large N c framework has not yet 
been reached. 



III. NP QCD VACUUM MODELS AND CASIMIR SCALING 

As it is mentioned in the Introduction, up to the authors knowledge there is no example of NP QCD vacuum 
model (with except of Gaussian stochastic model, see below) which does not violate CS in "naive" sense. 
All considered microscopic models need some kind of special adjustment in order to achieve approximate CS 
behavior of the static potential. We illustrate this statement taking as examples three popular pictures of 
the QCD vacuum - instanton liquid see ]2(| for review; condensate of abelian-projected monopoles ||27| . 
see ^8[^9) for review and center vortex scenario J3(J, |3l]j3^| . We will briefly discuss the confining string 
picture, D = 3 Georgi-Glashow model, MIT bag model and dimensional reduction scenario. 

To begin, let us make an important remark about the role of gauge invariance. As it was already mentioned 
in the Introduction, the property of CS is essentially nonabelian and therefore to correctly reproduce it a 
NP QCD vacuum model should take into account nonabelian degrees of freedom in a proper gauge-invariant 
way (roughly speaking, the sum in expressions like (||) must be taken over all N% — 1 generators). However, 
it is not sufficient since even manifest gauge-invariance of a model does not guarantee it will demonstrate 
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CS law. Nevertheless we find it instructive to separate the discussion of the models where some particular 
gauge-fixing procedure plays crucial role from those where it does not. We will see that pattern of CS 
violation is different in many respects for these two groups of models. 



A. Models with the preferred choice of gauge 



As the first example of a model, based on NP field configurations we mention here the model of dilute 
instanton ensemble [p5| , p6| . Strictly speaking this model is not expected to give reasonable predictions for the 
static NP potential since it lacks confinement and hence some essential part of NP physics. It is instructive 
however to investigate CS in this framework since it clearly illustrates the pattern which is common for all 
set of models based on ensemble of classical field configurations. 

The instanton ensemble is characterized by averaged instanton density n — N/V, where N is number 
of instantons and antiinstantons and V - 4-volume, and mean instanton radius p. Phenomenologically one 
chooses n « 1 Fm~ 4 and p « 0.35 Fm, in this case the "diluteness parameter" or packing fraction up 4 /N c is 
much less than unity, justifying the small density expansion in the model. 

The static potential in the dilute instanton gas was calculated in [B3| , see also p3] . In the approximation 



used in 33 34] each instanton contribution to Wilson loop is calculated independently, the potential at the 



leading order in density is given in this case by: 

oo 

V{R) = i J dp ^ J d s z Tr [1 - W(z, p, R)] (12) 
o 

where the instanton size distribution is such that / °° dn(p) — n ; pdn(p) = p n and the numbers 

of instantons and antiinstantons are assumed to be equal: Ni = Nj = N/2. The potential is normalised 
according to the condition V(R = 0) = 0. In case of static quarks belonging to higher representations eq. 
( |l2| ) should be generalised [M, with the following result: 

oo 

JeD 

and the function Fj(x) is given by some cumbersome double integral whose exact form can be found in j3jj ]. 
Here d(D) = D is the dimension of the representation D (not to be confused with the Casimir ratio do) and 
sum over J = 0, |, 1, .. goes over all SU{2) multiplets for decomposition of the given SUiZ) representation 

with the corresponding weig hts. One finds d{D) ■ C D = § £ J(J+ 1)(2J+ 1) and also £ (2J+1) = d(D). 

JeD JeD 
At small x the functions Fj(x) ~ x 2 , while at large x the functions Fj(x) tend to J-dependent constants 
p4| . Numerically one finds at small distances 

V(R) = 1.79 • 7 i? 2 -e,D+e>(i? 4 ) (14) 

where 7 = irpn and numerical coefficients 6d for D = 3, 8, 10 are given by 

e 3 : e 8 : do - 1 : 1-87 : 3.11 

instead of Casimir scaling results 1 : 2.25 : 4.5. Similar situation takes place for the large distance 



asymptotics of the instanton-induced potentail. It violates CS at the level of 20% . (seeJIO 



One-instanton approximation used in (12), ( |13| ) was recently examined on the lattice |35| in order to study 
the instanton density dependence as well as CS of the instanton-induced heavy quark potential. It was 
found that at the distance 1.2 Fm the deviation of the dilute gas formula ( |T2"| ) from the exact expression 
(H) is about 10% for fundamental charges and 8% for adjoint charges at the density n — 0.2 Fm~ 4 , while 
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it increases up to ~ 20% for fundamental charges and ~ 50% for adjoint charges at the phenomenological 
density n — 1 Fm' 4 . Moreover, linear density dependence of the potential extracted from the Wilson 
loop (|12| ) breaks down for the densities ~ 0.5 Fm -4 at the distance 1 Fm. The ratio of the potentials 
obtained according to (|l^) is found in j3f| at different distances and densities (to be compared with the 
Casimir ratio d 8 = C adj /C F = 9/4 = 2.25) and it goes as 2.12 at R = 0.4 Fm, n = 0.2 Fm" 4 ; 1.91 at 
R = 0.8 Fm, n = 0.6 Fm" 4 ; 1.49 at R = 1.2 Fm, n = 1.0 Fm" 4 . The ratio of the potentials extracted from 
(|J) has been found in |55| to be equal to 2 with an accuracy of a few percents for all explored distances. 

All these results make it clear that instanton ensemble cannot be a natural framework for explanation of 
CS. It also demonstrates that if one is considering field configurations on the lattice which by some reasons 
are mostly made of instantons (as it happens, for example, after the cooling procedure), CS does not hold 
in such vacuum. It would be interesting to observe directly on the lattice the CS breaking in this case. 

Let us mention another interesting problem in this respect - the pattern of CS in SUSY QCD. One should 
expect the role of instantons in vacuum dynamics of M = 2 and even M = 1 supersymmetric Yang-Mills 
theory to be much more transparent than in real world QCD. On general grounds one expects the simpler 
structure of perturbation theory and hence the absence of higher Casimirs in the perturbatively induced 
potential. On the other hand, coherent NP vacuum made mostly of instantons and antiinstantons should 
produce significant violations of CS for the NP part of the static potential. 

Another aspect is related to elementary fc-strings, which can be formed between static sources for the 
gauge groups SU(N C ) with N c > 3. Such strings are stable and characterized by the string tension crj. (see, 



e.g. P6J). Recent lattice studies in this direction [E3 24 37f| are of much interest. They are aimed to test CS 



in a different setup from that of j6j,|8| and compare the results with alternative predictions, e.g. with the 
so called MQCD conjecture |38|, see also |36|. This project is far from being complete and one needs more 
accurate data to make reliable conclusions, In particular, MQCD predicts <7fc cx sin(7rfc/iV c ), which gives for 
Ci/pfund the ratios 1 : 1.414.. : 1.618... for N c — 3,4,5, respectively. On the other hand CS hypothesis 
predicts that this ratio should scale as 1 : 4/3 : 3/2. It would be very interesting to test these predictions 
in further lattice calculations (see |E3 and references therein). 

Coming back to the real QCD vacuum, one can guess two possible explanations of the fact that instantons 
do not destroy CS. Either instantons are strongly suppressed in the real (hot) QCD vacuum (as it was 
observed in different respect in jjjj) while they are recovered by the cooling procedure. Or else instanton 
medium is dense and strongly differs from dilute instanton gas, in such a way that higher cumulant compo- 
nents of such collectivized instantons are suppressed. Interesting to note, that linear confinement missing in 
the dilute gas, would be recovered in this case. 

Now we come to another NP QCD vacuum model, based on the center dominance idea. An important role 



of the center of the gauge group was stressed by 'tHooft in 30 . A center vortices model tries to understand 
confinement in terms of interaction of the current forming the Wilson loop with the topologically nontrivial 
field configurations charged with respect to the center of the gauge group. A center vortex is a topologically 
stable field configuration with the topology of the surface, which carries (chromo)magnetic flux and interacts 
with the external (chromo)electric current in the following way: 

W(C) -> exp (2mnL(C, S) /N) W(C) (15) 

where n — 1,..,N — 1 and L(C,S) represents the linking number (in four dimensions) between Wilson 
contour C and vortex surface S. The idea of center dominance implies that dominant contribution to the 
string tension comes from the quantum fluctuations in the number of center vortices linking the loop. The 
model predicts the following static potential (the reader interested in the details of the center vortex model 
is referred to ]3l| , p2[ and references therein): 

N-l 

V D (R) = -J2 lo S (! - /n[l - ^ G D (c%{x))]) (16) 
n=i 

where the function Gd((^c( x )) ^ s defined as Gc(op(i)) = (l/d(D)) Tr exp[iaH] and {Hi} is the set of 
generators from the Cartan subalgebra. The function ot^{x) represents the corresponding solid angle and 



depends both on the Wilson contour C and on the position of the vortex center x. The parameter f n 
determines the probability for the vortex to cross with the contour. The expression ( filf) is derived in the so 
called fat center vortex model, where each vortex has finite thickness. This parameter is crucial for potential 
to be nonzero at intermediate distances in the case of zero triality representations: in the limit of vanishing 
thickness adjoint (and all zero iV c -ality representations) loop has no interactions with center vortices. So 
the introduction of thickness leaves a window for CS, which is essentially perimeter- type effect in this model 
p2| . The price to pay, however, is the contradiction between (|T(i| ) and (||), in particular, the potential in the 
adjoint representation is not linear (and even not a convex IffOf ) function of the distance, despite at some 
interval of distances it can be approximated by linear term and its flattening is going smoothly with the 
distance in this model, to be compared with (||). 

The exact value of the CS violation is strongly model dependent in this approach which makes it difficult 
to put stringent bounds on the parameters of fat vortices model based on CS. For example, in the SU(2) case 
the deviation of the center vortices induced potential from CS behaviour is about 30% for j > = 1 (adjoint) 
and ~ 80% for j = 3/2 at large distances |$2], but presumably it can be made much smaller by the proper 
adjustment of parameters. 

Lattice measurements of the potential for different representations of SU(3) were performed in center 
vortex model in J40] | . In is instructive to consider normalised Casimir ratio 

1 Vn(m 

io{R) ~ VAR) (17) 

In the theory with exact CS one expects £(i?) = 1. If we take for R = Rq some fixed point, for example 
such that RqVf{Rq) = 2.5 (this point corresponds to the edge of the region, where quasilinear asymptotics 
of Vf{R) begins), we get from the results (!(]] for some of the representations: 

&(Ro) = 0.98 , br(Ro) = 0.83 , £ 8 (2i? ) = 0.86 , Zi 5a (2R ) = 0.63 , £ W (3R ) = 0.47 

These numbers clearly demonstrate the level of the CS violation in the fat center vortex model. It was 
suggested in |}2| that the model can be tuned in order to explain properly the CS at intermediate distances, 
in particular, by adjusting the vortex profile function. At the same time, the initial root for CS violation lies 
in this case in the lack of the gauge-invariant formulation of the model. It would be interesting to understand 
in this respect the gauge-invariant contents if any of the notion of center dominance. 

It is of interest to discuss the situation which takes place in another popular confinement scenario - dual 
Meissner effect in the abelian-projected formulation of QCD |27],|53|]. We refer the reader to the review papers 
[^9|,^8|, where the details of this approach are extensively discussed. It is worth noting that at the present 
moment the dual Meissner confinement scenario is often considered as the most reliable microscopic picture 



of confinement and the method of abelian projections (in line with the abelian dominance hypothesis |54|) as 
the most adequate language for its description. As a warm-up example, let us consider first D — 3 Georgi- 
Glashow model. The SU (2) model contains nonabelian vector field with two components A^ acquiring 
nonzero mass in the Higgs phase, while the third component remains massless and provides dominant 
contribution to the Wilson loop. Existence of monopoles in this model leads to area law ]57j]. The CS 
property of Georgi-Glashow theory in comparison with abelian-projected QCD was discussed in Q. In the 
former case one has 



(Wj(C)) = (rr 3 Pcxp j ig J 4TM^J J ~ (tt 3 exp Up J A^dz, 

(exp ligm j A^zA ) (18) 



2j + l 

The average for large contours is dominated by m — or m = ±i for 2j even and odd, respectively. It is clear 
therefore that e.g. zero charge part of diagonal component of adjoint external current remains unconfincd 
and hence the model is CS violating. 
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In terms of gauge-invariant correlators given by (H^ (ki)..Hp n (k n )) where one defines abelian field strength 
Hp = e^pcf)" F a U p/mw , where 4> a is adjoint scalar Higgs field and mw - the mass of components of 
full nonabelian gauge field, the vacuum of Georgi-Glashow model made of 't Hooft-Polyakov monopole- 
antimonopole pairs is coherent and not stochastic. Indeed, nonperturbatively generated mass scale M ~ 
exp(— em^/2e 2 ) characterizing exponential falloff of the correlators in x-space enters string tension in highly 
nontrivial nonlinear way, while extracting confining part of the two-point correlator (Hp(k)H v (—k)) one 
gets simply D(k) — M 2 j(k 2 + M 2 ) Q. It is worth saying that nonperturbative physics of D = 3 Georgi 
Glashow model is investigated on the lattice in details (see, e.g. [^8| and references therein) and conclusions 
of theoretical analysis performed in J57j have been supported. 

We have learned from this example that Higgs mechanism may be a source of CS breaking since it makes 
some gauge degrees of freedom massive without destroying of gauge invariance, while other ones remain 
massless. This pattern is supposed to occur in QCD at finite density and we will briefly discuss it at 
the end of Section IV. Now we are coming back to abelian-projected QCD. Application of the method 
of abelian projections crucially depends on adopted gauge fixing procedure, so one could expect problems 
with the explanation of CS in this framework. Indeed it happens to be the case ]55| . Abelian dominance 
hypothesis |Q in QCD states, roughly speaking, infrared dominance of diagonal components A 3 ,^4 8 of the 
full nonabelian gauge field A® and hence suggests to omit nondiagonal parts in (|l|) in QCD just as it is 
done in (^). Doing so, one immediately obtaines deconfinement of adjoint charges at all distances in plain 
contradiction with the lattice data [j58j . It indicates that the dynamics of long-range QCD is not like that of 
duplicate compact QED (in contrast with Georgy- Glashow model, which behaves as compact QED at the 
considered limit) and color gauge symmetry is not broken by some composite Higgs field (as it effectively 
happens in abelian dominance scenario in QCD). 

There exist a few approaches in the literature aimed to cure this problem in the framework of abelian 
projection. All of them try to take into account the nonabelian degrees of freedom in more or less sophis- 
ticated manner. In particular, the reader is referred to |56| where the m = component of adjoint source 
is effectively dressed with a virtual cloud of charged degrees of freedom and hence is able to interact with 
neutral "photon" . Adjoint string tension arises in this picture at intermediate distances due to interaction of 
diagonal abelian projected gluons with the part of the adjoint source doubly charged with respect to the Car- 
tan subgroup. If one naively omits the corresponding Faddeev-Popov determinant it gives a a dj = 'l&fund, 
i.e. just the square of the abelian double charge. It is expected that loop expansion of the determinant 
produces terms, correcting the above behaviour to the Casimir scaling ratio, but it has never been shown 
explicitly. 

Another approach was discussed in [ p9| . The results of |5£j are based on the so called Weyl symmetric 
formulation of effective dual abelian Higgs model (DAHM) with the direct summation over S£/(3) root 
vectors. The following prescription for the external current interaction term was adopted (see also f34j) 

h^ a A a ^ - ^(T 3 Al + T 8 ^)^ (19) 

with subsequent recipe V>(:x;)7 M T 3 ' 8 -0(a;) — > Q 3 ' 8 • § c &{ z — x)dzp. The above ansatz is different from ( |l8| ) 
since it implicitely assumes some kind of averaging over different color components of external current and 
consequently produces (contrary to (|l8|)) area law for Wilson loop in all representations. It was found in 
this framework in [ [59| that CS favours different values of Ginzburg-Landau parameter for different represen- 
tations, lying in the Type II range k — ms/ma = 5^9. For k = 1 flux tubes become noninteracting and 
one has flux counting rule (JDl&fund = A* + v i where /i, v are Dynkin weights of the SU(3) representation 
(see |39| and references therein). 

The formalism proposed in |p0| deals with some lattice motivated reformulation of QCD, where a special 
kind of averaging over different choices of [£/(l)] subgroups of the original ST/ (3) is implemented instead of 
having one fixed choice. The situation with CS in this model if realistic values of Ginzburg-Landau parameter 
k are taken (e.g. from fitting of QCD string profile by DAHM induced formulas [j6l]| ) , is still unclear. We 
refer the reader to the review J|8| where this set of problems is discussed. 

Concluding the discussion of abelian projected theories let us mention that it could be of interest to 
calculate the potential in approaches advocated in [p6[j59|j60|| at all distances and for different representations 
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of SU(N) and directly compare it with the results of [[7|,|3|. As it is clear from the above discussion the 
problem to reproduce CS is still far from being closed for abelian projected theories. 

All the examples discussed above clearly demonstrates crucial importance to count properly all relevant 
nonabelian degrees of freedom. Any NP QCD vacuum picture analysed so far lacking manifest gauge- 
invariance had encountered internal difficulties explaining CS. In the next section we will consider another 
class of models where the nonabelian gauge invariance is kept intact from the beginning. 



B. Gauge-invariant models 



We start with MIT bag model |4l|] and refer the reader to papers 42 where the problem of the CS in 
this model was addressed. It was advocated the following relation between string tensions in the MIT bag 
model: a^jap = \fdr>. It is easy to see, that this relation is in contradiction with the data from at the 
level of a few dozen standard deviations. Qualitatively, it comes from the fact, that bag model misses the 
string, a crucial ingredient of modern picture of NP QCD, based on lattice results. The string in the static 
charge-ant icharge case has two ends, with the factor gTp attached to each one, and fixed, representation- 
independent radius (see below), so the total contribution to the energy in proportional to the charge squared, 
i.e. Cd- At the same time, in the MIT bag model the charge enters the energy linearly. 

Another family of models to be considered here are the so called "confining string" models. Attempts to 



build a theory of the confining string have begun more than 30 years ago (see e.g. 43 ). Till the present 
moment however there is no selfconsistcnt quantum string theory applicable to QCD string, those formation 
in the confinement phase is clearly seen on the lattice. Moreover, such theory is absent even for abelian 
Higgs model, where the strings exist as classical solutions (Abrikosov-Nielsen-Olesen strings) (see ^] and 
references therein) . As an example let us consider here the static potential induced by quantum Nambu-Goto 
string. 

The action of the model is proportional to area of a surface bounded by the static sources worldlincs. 
Quantum dynamics of this surface produces additional contribution to the confining potential besides the 
leading linear term [ [f~4"|| : 

a R^aR-^^ + ... (20) 

where the term — 7r/(12 • R) will be referred to as the String Vibration (SV) term. Despite the Nambu-Goto 
string model cannot be rigorously defined in D = 4 (see e.g. ji^j), and, in particular the expansion of the 
r.h.s. of (|2(]) is valid at large distances R only, it is instructive to look whether or not the lattice data support 
the existence of such term. It is also worth noting, that the dimensionless coefficient — 7r/12 is universal and 
determined by the only two factors: target space dimension and the chosen string model. Having both 
factors fixed, it cannot be freely adjusted. Assuming op, = dp>op, it is easy to see, that the Nambu-Goto 
SV term violates CS of the potential, as it was noticed already in ||. 

It is a nontrivial task to separate the contributions of the discussed sort in the confining potential as it 
is because these corrections are essentially large distance effect, where they are subleading. But they have 
to become pronounced in the expression (Q) due to scaling violation. Namely, assuming that SV term is the 
only one violating CS, one obtains 

V F (R) V F (R)Ro 12 R y ' 

where the dots denote terms, omitted in ( |20| ) and Rq is some arbitrary scale, for example Sommer scale. 
On general grounds one should expect, that string picture works at distances R l/y/ap since no other 
dimensionful parameter enters the Nambu-Goto string action (in particular, string in this model is infinitely 
thin). We expect therefore, that the data HQ] allow to extract the possible contribution from CS- violating 
SV term at intermediate distances ~ 1 Fm. The l.h.s. of ( |2~j| ) is compatible with zero at the level of la from 



11 



the data U while the r.h.s. of (|l|) is nonzero and is rising with do- However the statistical errors of the 
expression on the l.h.s. at such large distances are too big to make definite conclusion about ruling out the 
SV term of the discussed form. Additional accurate measurements are needed for this purpose. Notice, that 
the sign of the Nambu-Goto SV correction is opposite to what was found for CS violation in || ■ 

One way to combine CS and string dynamics was proposed already in || - just to multiply all the 
potential (pp| ) with dp, accounting for "Cjj elementary fluxes" attached to the static sources. It should be 
stressed that it implies different physical mechanisms responsible for creation of the string and its quantum 
fluctuations and presents actually a model, different from Nambu-Goto string. The situation with direct 
lattice measurements of SV corrections is not yet clear. While the authors of [0 claim the disagreement 
between hybrid spectrum and the string picture, but there are also evidences in favor of SV term [^6|. The 
question certainly deserves futher study. 

From theoretical point of view it stays to be proved that the QCD confining string and the simplest bosonic 
Nambu-Goto string model belong to the same universality class, see e.g. recent discussion in and the 
theoretical background of (|2(]) is not yet clear. There are also arguments against this possibility (see, e.g. 
J43"| ) . The theory of the QCD confining string - whatever it will be - must explain the observed Casimir 
scaling of the potential at intermediate distances. There are arguments [ |47| that the confining string in the 
strongly coupled continuous Yang-Mills theory posessess two limiting regimes: CS respecting, which occurs 
if the area of the surface bound by the Wilson loop is smaller than some critical area, and flux-counting 
regime, which starts to dominate for asymptotically large surfaces. It could happen that strings created 
by the sources in higher representations would break at smaller distances than those where the quantum 
fluctuations of the string become visible. In this case CS region will cover all the scales where the string 
exists as a stable object and no room to detect any significant contribution from SV terms will be left. This 
problem is absent for stable fc-strings on SU(N C > 3) (see above), and it is reasonable to expect CS violation 
at large distances in this case caused just by SV corrections, i.e. by quantum dynamics of the string. 

Let us mention two regimes of the SU(N) gauge theory where the CS property is under theoretical control. 
The first one is two-dimensional Yang-Mills theory. In 2c? case the theory has no propagating degrees of 
freedom and the static Coulomb potential which is linear in two dimensions is exactly proportional to 
Cd- Another example is the strong coupling regime, where the partition function of the Ad lattice Yang- 
Mills theory is dominated by the two-dimensional surfaces with the bare string tension proportional to the 
quadratic Casimir. It is worth reminding that the original motivation for invention of the CS hypothesis 
H was based on the dimensional reduction scenario which assumes that there exist nontrivial relation 
between NP observables, e.g. string tensions in 2d and Ad Yang-Mills theories, if proper identification of the 
couplings is done. QCD in 2d and Ad have indeed much in common while also some important differences. 
One of them to be mentioned here is the Lorentz structure of NP linear confining potential which is vector 
in 2d but scalar in Ad. We will briefly discuss the relation of this picture with Gaussian dominance in Section 
V. 

In the gauge-invariant NP background field formalism (see review [||) the CS property has two possible 
interpretations. To understand the corresponding physical pictures, let us formulate some basic ideas of field 
correlator method. The approach decomposes the gluon field into perturbative field and nonperturbative 
background with the latter to be taken into account by means of gauge-invariant irreducible correlators (cu- 
mulants) ((Tr <I>(xo, xi)F{x\)<&(xi, xq)...<&{xq, x n )F(x n )$(x n , xo)}}- In particular, the Wilson loop average 
in the " 34" plane admits following expansion (see for details) 

(W D (C)) = (rmPexp L J A a ^T a dz)j ^ = (ri-oVexp L J F a (u)T a da(u) 

OO n OO 

= Tt d T x exp / i n {{F{u^)..F{u^)))da{u {l) )...d(j{u {n) ) = Tr D exp ^™A^ } [S] (22) 

n=2g n=2 

The nonabelian Stokes theorem j4S| and the cluster expansion theorem have been used in deriving (|2 
Here F(u)da(u) — ^(xQ,u)E^(u)T a ^(u,xo)da34 : (u), where $ is a phase factor, 
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$(x ,w (fc) ) = Pexp 



(23) 



«W and xo are the points on the surface S bound by the contour C = dS. The double brackets ((...)) denote 
irreducible Green's functions proportional to the unit matrix in the colour space (and therefore only spacial 
ordering V x enters (22)). Since ( 24) is gauge-invariant, it is convenient to make use of generalized contour 
gauge which is defined by the condition $>(xo, u^) = 1. 

In the confinement phase vacuum is disordered in a sense that averages from the r.h.s. of (^) develop 
finite correlation length T gi which for the lowest two-point correlator was found on the lattice to be rather 
small T g « 0.2 Fm in the quenched approximation for SU(S) (to be compared with the old stochastic 
proposal of pOfl). The Gaussian dominance Q implies that the dominant contribution to the potential (and 
most of other observables) comes from the lowest two-point correlator (FF) if the integration surface in ( p2| ) 
is taken as minimal. It is easy to see that 

(Tr D F(1)F(2)) = (F«(1)F«(2)) = ^ (F»(1)F°(2)>, (24) 



so Gaussian approximation provides exact CS. Notice that (24) holds true since parallel transported field 
strength tensors F(u) — $(xo, u)F a (u)T a <fr(u, xq) from fl22] ) obey the same standard commutation relations 
as bare F^(u)T a do. 

The CS property of Gaussian vacuum depends neither on the actual coordinate profile of the potential 
nor on that of the correlator, in particular, it holds for any value of T g . It is true also for representations 
where the linear potential is just some kind of intermediate distance behavior and changes profile at larger 
R. The coordinate dependence of the potential, not directly related to CS, can be analysed at the distances 
small enough not to be affected by screening effects. 

The average (Wd{C)} does not depend on the choice of the surface S in (p2|). At the same time the 
particular choice of S determines the weight of integral contribution of each cumulant A^ [S] to the Wilson 
loop average, and this weight is obviously S'-dependent. However this dependence is cancelled out if the 
whole sum in the r.h.s. of ( p2| ) is taken pS| ]. It is clear therefore that for surfaces strongly different from 
the minimal one there is some kind of cancellations between higher terms, while for the minimal surface 
CS might be either a result of cancellations between higher Casimir contributions from higher cumulants 
or just follow from the Gaussian dominance (terms with n > 2 in ( p2| ) give small contribution). We will 
discuss below possible ways to check that. In either case, we should stress again that Gaussian dominance 
is understood throughout this paper in the following sense: if S — S m i n is the minimal surface, A^, [S„] 
in (|2|) is considerably larger than all other terms. The gauge-invariance of A^[S] makes the notion of 
Gaussian dominance gauge-invariant as well. 

IV. WHAT CAN WE LEARN FROM THE CS? 

In this section we discuss other consequences of CS and try to argue that some further detalisation of the 
CS property in lattice calculations might be of much use for our understanding of the strong interaction 
physics and the mechanism of confinement. 

A. Three-point correlators 

Only two-body potentials of the meson type have been discussed till the present moment. It is also of 
interest to study three-body system, which in fundamental and adjoint cases corresponds to baryons and 
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tree-gluon ghreballs, respectively. The tree-body potential for fundamental sources V 3Q (R) is denned by 
expression (g) where W^ 3 ^ is taken to be 

WW = e ct/37 $^(C 1 )^(C 2 )$^(C 3 )e' ,CT " (25) 

and the contours Ci, 62,63 are formed by static quark trajectories. Indices of the antisymmetric tensors 
e Q /3 7 run from 1 to 3. For ajoint sources the corresponding expression is given by 

w^ G) = n abc $ ad (d)$ foe (c 2 )$ C5 (c 3 )r! deg (26) 

where the tensors f2 a h c determine the colour structure of the 3G state corresponding to the minimal energy 
and adjoint indices run from 1 to 8. Let us take spacial configuration of the system to be equilaterial triangle 
with the side R. Naively one might assume the same CS law for the 3-body case, i.e. V 3G (R) /V 3 ® (R) = 
Cadj I Cf ■ 

However, one has to conclude that the CS hypothesis in 3-body case should fail. The reason lies in different 
geometry of the 2-body and 3-body systems. Indeed, in case of two static sources the minimal area or the 
surface bound by the Wilson contour is provided by the plane, one and the same for all representations. In 
3-body case in adjoint (and higher) representation the string may go either along Mercedes-star configuration 
(with the formation of a string junction, Fig. la) or as triangle configuration (see Fig. lb) and one can easily 
argue that the energy of the latter variant is lower (notice at the same time that there is no such choice in the 
case of fundamental sources). Indeed, for the static (i.e. noninteracting and nonvibrating) strings connecting 
adjoint sources one has V^ } (R) = ap • 3 R for the triangle configuration and V 3G (R) — a a dj • 3 • (i?/v3) = 
(9/4)<r F • V3R > Vl G (R) for the Mercedes-star case, where CS relation a a dj = (9/4)ctf was used. For higher 
representations one can imagine more complicated geometry of strings, for example with more than one 
junction. Therefore we do not expect that lattice analysis of the 3-body potentials for higher representations 
will indicate CS in the same fashion as it does for 2-body potentials. In particular, it is interesting to check 
the ratio V 3G (R) /V 3G: (R) ~ V3 which would indicate different geometry of the strings (and at the same 
time indirectly support CS law) in the two cases. 



B. Stochastic vs coherent 



It was shown in the previous section how the gauge-invariant picture of Gaussian dominance explains CS 
of static potential. It was also stressed that the CS property on the other hand does not imply Gaussian 
dominance since one can imagine an alternative scenario (to be referred to as 'fine-tuning' picture) where 
all higher terms in the expansion (^2|) strongly contribute to the lowest Co-proportional term together with 
Gaussian correlator while all contributions with higher powers of Co and higher Casimirs effectively cancel 
each other. Therefore despite CS property seems to be a very strong evidence for Gaussian dominance, 
the latter stays to be checked explicitly. See |54| in this respect where the quartic field correlator was 
measured and evidences for its smallness were found. We propose here several types of field correlators 
which distinguish coherent and stochastic (i.e. Gaussian) pictures. The first one is defined as follows: 

(Tr $(z 0) x)F lxv {x)^{x, xoMx , y)F pa {y)<S>(y, x )) = D$ pa {x -x ,y- x ) (27) 

The phase factors are defined according to ( p3|) with the straight lines connecting edge points as integration 
contours (see Fig. 2a). The xo-dependence of ( |27j ) is related to the contribution of correlators higher than 
Gaussian, but respecting CS. Gaussian dominance implies smallness of this contribution. Notice, at the 
same time that ([27|) is exactly proportional to the quadratic Casimir by construction. We take all fields in 
the expression (|27|) in the fundamental representation for the sake of simplicity. Using standard relations 
for the phase factors (see, for example, pSy ) one gets: 

l 

d f 

D $pcr = i 9 dss ( x ~ £o)/3(Tr F piV {x)[F [ja {x + s(x - x))F pr7 (y)}) 

o 
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ig / dss(y-x )p('TrF liV (x)[Fp a (y + s(x -y))Fp <T (y)]) 



(28) 



with the explicit contributions from the 3-point correlator to the r.h.s. of ( f28|) (we omit <£>'s for brevity). 
One can also define two typical lengths, characterizing the behavior of ([27]) in different regions of parameter 
space. It is physically clear that the situation when \x — y\ ~ \x — Xq\ + \y — xq\ should be distinguished from 
the case with the distant reference point \x — y\ <C \x — xq\ + \y — Xq\. The longitudinal correlation length 
in corresponding to the variation of x is to be compared with the transverse one L±, which parametrize 
the change of ( f27| ) with respect to xq in both regions in lattice simulations. It is also possible to choose 
rectangular geometry, it corresponds to the correlator 



(Tr F M „(0, 0)$(0, 0; z, 0)<!>(z, 0; z, T)$(z, T; 0, T)F pa (0, T)$(0, T; z, T)$(z, T; z, 0)$(z, 0; 0, 0)) 



(29) 



where all integration contours are again straight lines (see Fig. 2b). The z-dependence of (|2^) is to be 
explored. The correlators of this form enter different expressions when one quark is heavy and can be 
considered as a static source. In Gaussian dominance scenario one takes into account T-dependence and 
neglects z-dependence of ((H) |6ql . 

The second type of correlators directly distinguishing stochastic and coherent pictures is given by multi- 
point irreducible field strength correlators. The simplest nontrivial one is that of the fourth order (see Fig. 
2c, where the choice of contours for four-point correlator with all arguments on straight line is depicted). 
One has to consider the following expression 



xi x 2 x 3 

&d ] [S}= J da{x x ) J da{x 2 ) J da(x 3 ) J da{x 4 ){Tr D F(xi)F(x 2 )F(x 3 



)F(x A )) 



■ ri 

da( Xl ) J da(x 2 )(Ti D F(x 1 )F(x 2 ))- 



(30) 



where we omit Lorentz indices and phase factors $>(xi,Xj). Notice that in confinement phase where field 
correlators develop finite correlation length, each term in the r.h.s. of ( ^0[ ) is proportional to R 2 at large 
distances, while their difference and hence the l.h.s. of (|3(]) is linear in R. 

The ordered integration goes over the surface S, parametrized by x(s,t) : x(0,t) = x ;x(l,t) G dS 
according to 



da^{x 2 ) 



= / dsn 



dtn 



( dx^{s 2l t 2 ) dx^(s 2 ,t 2 ) _ dx^(s 2 ,t 2 ) dx%(s 2 ,t 2 ) 
\ ds 2 dt 2 ds 2 dt 2 



It should be clear from the previous discussion that in Gaussian scenario such surface-dependent contribution 
A^' [S] is assumed to be small in comparison with the two-point one 



da{x x ) / da(x 2 )(Ti D F( Xl )F(x 2 )) 



(31) 



if the surface S is the minimal surface. There exist lattice indications |64| that this is indeed the case. It 
would be of considerable interest to explicitly estimate ( |30|) in comparison with (^l|) on the lattice using 
powerful present day lattice technique and also to study the CS behaviour of (|30|). 

The suggested probes being particular examples of path-dependent gauge-invariant field correlators are 
more subtle than Wilson loop averages and can be effectively used to discriminate coherent and stochastic 
pictures in direct way. 
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C. NP effects at small distances 



It was recently suggested on both theoretical |6(| and lattice p7| grounds that leading NP contribution to 
the static potential at small distances linearly depended on R, instead of the quadratic dependence naively 
suggested by conventional OPE: 

V D (R) = -C D ^+a D -R (32) 

where the dimension-two coefficient &d should be distinguished from the asymptotic string tension an in (g) . 
From the background perturbation theory point of view described in the previous section, this linear term 
comes from the interference between perturbative and nonperturbative interactions. In theoretical picture 
developed in |3(| the phenomenon was associated with formation of the so called infinitely thin string between 
charges. One can address the issue of CS in this approach. It is of interest to understand theoretically and 
on the lattice in this framework the relation between "thick confining strings" which provide confinement 
and infinitely "thin short strings" violating conventional OPE and presumably contributing to the linear 
potential at small distances in the light of CS in both the small and the large distance regions, found in 0. 

There are arguments to be mentioned |6^] that the quasilinear behavior of the potential at small distances 
observed in [p7[ could be a manifestation of the phenomenon of freezing of perturbative coupling at large 
distances. It is interesting therefore to check CS in this setup directly, namely to extract not only ap, as it 
was done in (67) but also the same quantity for higher D in analogous way. The linear scaling of &d with 
Cd would be in line with the proposals of |]66|| and of course with Gaussian dominance. Hopefully further 
lattice investigations will shed light on that. 



D. Interaction of Wilson loops 

Another lattice setup where the CS phenomenon can be successfully checked is the Wilson loop correlators. 
In the present subsection we start from the disconnected correlators of the form ([\ WoiiCi)} and come to 

i 

the discussion of the connected ones (Tr d II < J > u(Ci)) after that. The case of interaction between loops in 

i 

different representations is of particular interest. We consider the two-loop correlator in the regime where 
the sizes of both loops are much smaller than the distance R between them. One gets (see Fig. 3a) 

(W Dl (Ci) ■ W D2 (C 2 )) - (WnACr)) ■ (W D2 (C 2 )) ~ [C Dl ■ C D . 2 + ...] • cx ^^ (33) 



where M g i does not depend on the choice of representations D\,D2- The terms denoted by dots in ( |33j ) are 
proportinal to the higher powers of Cd; and should be suppressed if Gaussian dominance holds true. Let 
us mention the results from |69f| where some indications were presented that (^3|) holds true (for adjoint- 
fundamental correlator). See also |Q where the behavior (^) was obtained in the Gaussian stochastic 
model framework. 

Let us consider now gauge-invariant connected correlators. Lattice simulations demonstrate the phe- 
nomenon of confining string formation. The important problem of confinement dynamics is to study the 
string profile. This question can be addressed for the string attached to fundamental sources and to sources 
in higher representations as well. One's first instinct is to say that since for higher representations the string 
tension is larger the same should in some sense be true for the geometrical characteristics of the string, e.g. 
its radius. It was already discussed in the previous section that it is in fact not true. To see this, let us 
consider connected probe of the following form: 

d = (Tv^ D (C)^(x)^) _ 
*V (Tr*£>(C)) ( ' 
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where phase factors $^ connect the plaquette P^(x) with the loop $d(C) (compare with (| 
Since in the continuum limit {P^u{x))qq — > a 2 {F^, v {x))qq where a is lattice link one gets 



c 



D 



S(C) 



+ 0(a 4 ) 



(35) 



where S(C) is the minimal area surface inside contour C . Omitted terms in the r.h.s. of (|35|) denoted by 
dots contain higher powers of Cd- It is seen at the same time that the string radius (d\ogpu/dx) 1 is 
representation-independent in Gaussian stochastic model (and, generally, for arbitrary CS respecting field 
ensemble), if the terms proportional to Co dominate in (pa). On the contrary in the absence of CS for 
( pM| ) the representation dependence of log needs not to be factorizable and radius of the string for such 
coherent ensemble would be representation-dependent. 

It is therefore gauge-invariant analysis of the QCD confining string profile for charges in higher represen- 
tations which could indicate both CS and Gaussian dominance in a way different from just static potential 
measurements. Notice, that in [[39 70 another (disconnected) type of correlator was used and the results for 
(|35|) should not be directly compared to that of those papers. 



E. CS for nonstatic potentials and hadron spectra 

Till the present moment we have considered only static potentials. From hadron phenomenology point of 
view, however, not static but dynamical potentials, which take into account the effects of quark motion are 
of interest. What is the meaning of CS in this case? Let us consider for simplicity gauge-invariant Green's 
function for heavy-light spinless system in the gluon field (see, e.g. JiJ) 



G(x, y) — Af ds / Vz exp 



-m 2 s 



da- 



1W 



Tr £>Pexp 



+i / Al(z,z 4 )T a d ZfJ 



P exp 



iji 



i / Al(6,t)T a dt 



(36) 



There are several different regimes of (36) depending on relative value of the light particle mass m and 
properties of the gluon ensemble, encoded in the weight used in the averaging procedure (..). In case of large 
to when the nonrelativistic expansion is valid, one has to distinguish two situations: the potential regime and 
the sum-rule limit. The energy levels get a correction in the latter case which in Gaussian stochastic model 
is proportional to Cd ■ Contributions of higher condensates are suppressed in this case by the inverse powers 
of to. In the potential regime the situation becomes less transparent and one is not to expect any exact 
scaling law with respect to Cd for the spectrum of the system, despite the potential obeys CS. Spin degrees 
of freedom provide additional complications. In particular, difference in spin structures precludes to identify 
the spectrum of the heavy-light system where spinor light particle carries fundamental color charge (which 
models D and B-mesons) with heavy-light spectrum where vector-like dynamical constituent carries adjoint 
colour charge (such system might mimic hypothetical hybrid excitations of heavy mesons). Indeed, the only 
direct fingerprint of CS for the real particles spectrum is presumably the fact that the lowest glueballs are 
believed to be heavier than the lowest mesons (we do not take chiral effects into account). Let us consider 
"hadrons" as bound states of massless and spinless quarks, without perturbative interactions. In this case 
the mass of a hadron is proportional to the square root of od - the only dimensionful parameter of the 
problem: M n = c n ■ ^od ~ VCd- Hence one might predict that the equivalent states in glueball and meson 
spectrum have ratio of masses M" /M 9< ? = \/9/4 = 3/2. Now include spins and perturbative dynamics. 
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Taking into account of Coulomb force changes this ratio, also spin-dependent forces are different since spin 
of gluon is twice that of quark, and finally quarks have (negative) self-energy corrections (of the order -0.25 
GeV per particle J7l| ) while for gluons selfenergies are forbidden by gauge invariance. All that leads to the 
masses calculated in potential model in f72|| which are in good agreement with lattice and experimental data, 
namely for center-of-mass of 2++ and 0++ glueball one finds M c 9 o 9 i L=0) = ±(M 9S (2++) + M^(0++)) w 2 

Gev, while m p = 0.77 GeV and similarly for M^ =1) = 2.66 GeV, while for mesons M^ =1) = 1.2 GeV. 
Ratio of glueball to meson masses in both cases is around 2.2 — 2.5, which is larger than 3/2 = 1.5. The 
main reason is quark selfenergies and different role of Coulomb forces: while for qq those are significant, for 
gg system as argued in J72| perturbative gluon exchanges are suppressed due to higher loop effects. 

F. CS and QCD sum rules. 

The OPE series which is the starting point of QCD sum rules contains two types of higher twist term: 
those which contain color-irreducible combinations of field operators, like (D k F^D n F a \) or ("ipa^F^t/j) 
and reducible operators like ('tptp'tptp) or (FFFF). For the latter the conventional strategy |73| is the so 
called vacuum insertions, replacing e.g. (FFFF) — > (FF)(FF). This replacement procedure was however 
criticized from the point of view of the instanton gas/liquid model, which does not suport good accuracy of 
vacuum insertions. Now after measurements done in [@j8| and subsequent analysis in [jTo| , p^| one can say more 
about accuracy of vacuum insertions at least in purely gluonic operators. Namely, such vacuum insertion 
is equivalent to neglection of higher irreducible correlators (cumulants), e.g. ((FFFF)) ~ (FFFF) — 
(FF)(FF) ps and since quartic and higher order correlators violate CS, their admixture is limited by the 
measured accuracy of CS, i.e. of a few percents. The same CS also shows (see Section III of this paper) 
that instanton model of vacuum fails to reproduce the correct behavior of static potential and hence the 
procedure of vacuum insertions for gluonic operators is fully justified with shown above accuracy. Roughly 
speaking, CS favours vacuum insertions. 

Another facet of that is also worth mentioning. From the spectral decomposition point of view the vacuum 
insertion procedure is based on the fact that there is a mass gap in the theory and a set of massive (colorless) 
states gives small contribution to the spectral decomposition of color-reducible correlator. On the other hand, 
vacuum is made of colorless weakly interacting dipoles in Gaussian picture and confinement (and hence the 
mass gap) is associated with the fact that these dipoles have finite (and actally small) size. The size of these 
dipoles is inversly proportional to the lowest gluelump mass (see above). On the contrary, the formal OPE 
limit is achieved by taking this size as infinite, T g — > oo in terms of FCM, when field correlators become just 
x-independent condensates. It is of interest to clarify this problem selfconsistcntly and to establish the form 
of the operator product expansion in the physical Gaussian vacuum, characterized by the finite correlation 
length. This work is in progress now. 

G. CSatT/0 and y, / 

Lattice simulations of QCD at nonzero temperature could provide additional information about the dis- 
cussed question. One should distinguish the Polyakov loop, i.e. trace of the phase factor taken along 
temporal axis from to (3 — 1/T, and usual Wilson loop, which is made of spacial links. It is known (see 
e.g. ^ and references therein) that at the point of deconfinement phase transition the magnetic components 
of the chromoficlds stay approximately intact while the electric components responsible for the nonzero 
string tension of the temporal Wilson loop vanish. In accordance to that we expect that the CS law for 
the potential extracted from the spacial loops not to be spoiled at T > and even at T > T c . As for the 
Polyakov loops, the CS has no simple meaning for the averaged traces of temporal links. At T < T c the 
Polyakov loops for nonzero triality charges vanish while the screening happens for other representations. 
Therefore there is no room for CS region in this case. After the deconfinement transition one comes back to 
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the dynamical CS and has to study the spectrum of heavy-light lumps with the heavy source taken to belong 
to the given representation. As it was argued in the previous subsection, such spectra need not necessarily 
obey CS. At the same time it is important to notice that lattice shows [ |74| one and the same deconfinement 
temperature extracted from Polyakov loops in different representations. This fact seems to have no direct 
relation with CS, at the same time, it is very natural in Gaussian dominance picture - it just indicates that 
the electric part of the two-point field strength correlator, responsible for nonzero string tension, vanishes 
above the critical temperature. The actual value of the deconfinement temperature does not depend on the 
properties of external current. It would however be very interesting if lattice simulations will indicate that 
CS is temperature-dependent, studying either spacial or temporal loops. In Gaussian dominance picture it 
will signalize that the dynamics governed by higher correlators is different at different temperatures. 

Another closely related problem is the determination of higher representation Polyakov loop behaviour 
when T — ► T c from above (see, e.g. [fT5[ , (7f|). The difference in critical exponents for higher representation 
Polyakov loops near the point of the phase transition presumably can be interpreted in terms of some 
bounds on higher cumulant contributions in (|2^). This set of questions in connection with CS will be 
discussed elsewhere. 

The situation at /i 7^ is much more subtle. It is supposed J77| that at high density the phenomenon of 
color superconductivity takes place. One should distinguish two and three flavor cases. In the latter case 
the so called color- flavor locking leads (see review |78 and references therein) to nonzero and equal masses 
for all eight adjoint colors of gluons, while in the two-flavor case vacuum state breaks global color symmetry 
SU(3) — > SU(2) in such a way that three gluons corresponding to unbroken gauge subgroup SU(2) remain 
massless while five others take nonzero and unequal (in low temperature phase) masses. The picture is 
therefore in some correspondence with the Georgi-Glashow model discussed above and it is reasonable to 
expect that CS is destroyed in the region of densities corresponding to two-flavor color superconductivity 
(and restored back again at densities where the third flavor comes into play) . Despite external Wilson loop 
seems to be rather artificial probe at high density and lattice simulations of N c — 3 QCD at // 7^ are 
unfortunately impossible at the moment, it is tempting to speculate that two-flavor color superconductor if 
exists might occupy the only region on the QCD phase diagram where the CS property is broken in explicit 
way. 



V. CONCLUSION AND OUTLOOK 

Accurate lattice data on the static potential for the sources in higher SU(3) representations strongly 
support the Casimir scaling (CS) hypothesis Most of the popular models of confinement encounter 
internal difficulties in attempts to explain this fact. We have shown in the present article that CS strongly 
supports the Gaussian dominance of QCD vacuum. Let us summarize the paths which we propose for further 
studies of this complicated dynamical problem. 

• To study CS of perturbative series beyond two loops both in continuum limit and on the lattice. 

• To investigate CS in the course of the cooling procedure. 

• To establish exact results for the static potentials between charges in higher representations in SUSY 
QCD and compare them with CS-respecting results of ordinary QCD. 

• Convincing evidences for the string breaking phenomenon will allow to establish quantitatively the 
bounds of the CS region. 

• CS might break at large enough distances were the quantum dynamics of the confining string starts 
to be detectable. Needless to say that it would be very interesting to observe the regime of vibrating 
confining string, if it exists. 
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• The set of lattice measurements which could significantly improve our understanding of CS and Gaus- 
sian dominance includes static 3-point correlators, nonperturbative potential at small distances, string 
profile, Wilson loop correlators in higher representations etc. 

• It should be stressed again that despite Gaussian dominance straightforwardly implies CS, the opposite 
strictly speaking is not true. By direct measurements of path-dependent field strength correlators one 
can check explicitly the accuracy of Gaussian dominance (i.e. comparing ( |30| ) and (j31~l)). 

Whatever the dynamical reason of CS is, its persistence strongly advocates that QCD vacuum is stochastic 
rather than coherent. It provides also some limits on the use of instanton gas/liquid model and other models 
of confinement based on classical solutions. As was explained above stochastic picture of QCD vacuum 
incorporates some features of the dual Meissner effect. 

In terms of field correlator method Gaussian dominance implies the dominant NP contribution to the 
gauge-invariant observables from the NP parts of the two-point field strength correlator. Physically it 
corresponds to the picture of the vacuum, made of relatively small white dipoles. The size of these dipoles 
which is natural to associate with T g controls the stochasticity of the vacuum, i.e. the spacial size of the 
domains where fields are coherent in the spirit of old ideas [Eoj. We need however more than that - indeed, 
Gaussian dominance implies the approximate ideality of this "dipole gas", i.e. relatively weak interaction 
between dipoles and as a consequence a small fraction of the higher "multipoles" in the vacuum. One can 
possibly understand this structure solving nonlinear equations for field correlators generalized by inclusion 
of perturbative and higher correlator terms. 

The discussed white "dipoles" live on the surfaces bound by Wilson contours and the relative integral 
weight of higher irreducible cumulants (" multipoles" ) depends on the shape of these surfaces despite total 
answer is of course surface-independent and the dynamics of the dipole ensemble is topological in this sense. 
Gaussian dominance means that one can find such surface S - the minimal one - where the dipoles effectively 
become quasifree. Notice the difference from the dimensional reduction scenario where the effective 2d surface 
is populated by colored vector particles - gluons - with no any interaction between them for arbitrary 
geometry of the surface. 

Finally, the greatest mystery which probably explains all other features is the small value of gluon corre- 
lation length T g (or, equivalently, large value of the lowest gluelump mass, which is around 1.5 Gev) in units 
of string tension. Taken together with even higher values for multigluon gluelumps, yielding the correlation 
lengths of higher cumulants, this might explain both the Gaussian dominance and the high accuracy of CS. 
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